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ABSTRACT 

An  effective  data  interface  between  a  Victoreen  PIP-400 
pulse-height  analyzer  and  a  Hewlett-Packard  9810A  calculator 
was  designed,  built,  and  tested.  A  calculator  progrcim  was 
written  which  enabled  a  research  group  studying  laser 
propagation  in  the  marine  boundary  layer  to  conduct  rapid, 
local  processing  of  scintillation  and  extinction  data. 
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I. 


INTRODUCTION 


The  potential  utility  of  laser  technology  in  the  areas 
of  line-of-sight  coinmiinication,  target  tracking,  and  wea¬ 
ponry  makes  an  understanding  of  the  propagation  of  laser 
beams  near  the  air-ocean  interface  imperative.  The  effects 
of  a  turbulent  atmosphere  on  laser  transmission  were  theor¬ 
etically  modelled  by  Tatarski  [Ref.  1]  et.  al.  and,  subse¬ 
quently,  extensive  field  research  in  the  terrestrial  boundary 
layer  has  substantiated  the  bulk  of  the  theory.  Research 
data  in  the  marine  boundary  layer  has  been  scarce  and 
inconclusive . 

Atmospheric  effects  on  laser  propagation  may  be  divided 
between  laser- induced  phenomena  and  those  caused  by  meteoro¬ 
logical  characteristics  independent  of  the  laser's  presence. 
The  first  classification  applies  only  to  lasers  of  sufficient 
power  to  heat  the  atmosphere  (thermal  blooming,  upwind 
drift) .  The  second  is  a  universal  grouping  affecting  the 
whole  electromagnetic  spectrum  regardless  of  power.  The 
atmospheric  effects  include  beam  refraction,  beam  spread 
(or  degradation  of  resolution) ,  absorption,  and  scattering. 
Refraction  may  be  further  categorized  according  to  the 
dimensions  of  the  atmospheric  whirls.  Large-scale  vortices 
cause  path  deviation  (beam  wander)  while  small-scale  vortices 
cause  time- intensity  variations  (scintillation)  and  degrada¬ 
tion  of  resolution. 
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The  Naval  Postgraduate  School,  Monterey,  has  been  involved 
over  the  last  four  years  in  a  continuing  effort  to  establish 
the  data  base  necessary  to  confirm  or  modify  existing  theory 
as  applied  to  the  air-ocean  interface.  Two  research  groups 
have  operated  simultaneously  but  independently;  the  optical 
research  group  and  the  atmospheric  turbulence  group  have 
attempted  to  correlate  determinations  of  (the  refractive 
index  structure  constant)  through  the  respective  analysis 
of  light  intensity  distributions  and  meteorological  data. 

The  optical  research  group  has  concurrently  investigated 
scintillation,  extinction  (scattering  and  absorption)  and 
the  atmospheric  modulation  transfer  function  (a  measure 
of  beam  spread  and  wander) .  Theses  on  the  effects  of  aperture 
averaging,  determination  of  an  atmospheric  extinction  con¬ 
stant,  and  determination  of  an  atmospheric  MTF  for  a  black- 
body  source  were  completed  over  the  same  period  as  the 
author's  [Refs.  10,  11,  12]. 

A  major  problem  for  the  optical  group  was  the  inordinate 
amount  of  time  necessary  to  process  scintillation  and  extinc¬ 
tion  data.  The  methods  employed  depended  heavily  on  operator 
analysis  of  raw  data  or  extensive  manual  input  to  computers. 
While  the  processing  theory  being  applied  was  appropriate, 
the  time  penalty  and  the  inherent  high  probability  of  error 
had  to  be  minimized.  The  author  concentrated  his  efforts 
in  this  area  employing  only  the  basic  results  of  laser 
propagation  theory . 
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II.  PROJECT  BACKGROUND 


Various  members  of  the  Naval  Postgraduate  School's 
faculty  have  been  conducting  research  on  atmospheric  pro¬ 
pagation  phenomena  since  1972.  The  first  experiments  inves¬ 
tigated  the  atmospheric  MTF  and  scintillation  in  the  basement 
corridor  of  Spanagel  Hall.  In  the  scintillation  experiment 
a  laser  beam  reflected  from  a  flat  mirror  placed  at  the 
opposite  end  of  the  137  meter  passageway  and  entered  a  pin 
diode  detector  located  near  the  laser  source.  The  detector 
output,  proportional  to  the  incident  beam  intensity,  was 
passed  through  a  log-amplifier  and  recorded  on  magnetic  tape. 
The  data  was  later  periodically  sampled,  digitized,  and 
accumulated  by  a  pulse-height  analyzer;  the  resultant  dis¬ 
tribution  was  Gaussian,  as  theory  had  predicted.  Refractive 
index  structure  constants  were  calculated  from  photographs 
of  the  distributions  by  measuring  the  curves ’  full  width  at 
half  maximum  height.  The  MTF  measurement  used  a  six  inch 
parabolic  mirror  mounted  as  a  Newtonian  telescope.  The 
image  from  a  laser  at  the  other  end  of  the  passageway  was 
scanned  across  a  slit  using  a  cam-driven  rocking  diagonal 
mirror.  The  output  of  a  photocell  was  fed  to  a  boxcar  inte¬ 
grator  to  generate  time-averaged  data,  which  was  then  recorded. 
Turbulence  in  the  hall  was  somewhat  controllable  and  detec¬ 
table  reduction  of  MTF  with  turbulence  was  observed. 
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In  late  1972  a  group  of  interested  researchers  discussed 
possible  environmental  experiments.  The  first  such  experi¬ 
ment  was  completed  10  January  1973  by  faculty,  staff,  and 
students  and  involved  shining  a  low-power  (mW)  red  laser 

O 

(6328  A)  across  a  4.3  kilometer  range  at  the  southern  end 
of  Monterey  Bay.  The  research  vessel  Acania  was  anchored 
at  the  approximate  midpoint  of  the  beam  path  to  take  atmos¬ 
pheric  turbulence  measurements  with  onboard  sensors.  In  a 
second  experiment  at  San  Nicholas  Island  (3  March  1973)  low 
power  laser  beams  were  transmitted  from  the  northwest  tip 
of  the  island  to  the  Acania  anchored  offshore.  The  qualified 
success  of  these  experiments  completed  with  equipment  already 
owned  by  the  school  was  encouraging  enough  to  prompt  the 
sxibmission  of  a  research  proposal  to  the  Naval  Ordnance 
Laboratory  (White  Oak,  Md.)  for  fiscal  year  1974. 

Experimental  apparatus  has  changed  significantly  since 
the  Navy  lab  started  sponsoring  the  research.  A  Cassegrainian 
telescope  system  with  an  eighteen  inch  parabolic  mirror  and 
an  elastically  mounted  scanning  mirror  (driven  electromag- 
netically  by  a  shaped  waveform)  was  used  in  the  most  recent 
MTF  measurements .  An  enclosed  trailer  was  built  to  transport 
and  protect  the  telescope  and  associated  equipment;  aluminum 
support  legs  that  extended  to  rest  directly  on  the  ground 
negated  the  previously  persistent  problem  of  telescope  motion 
due  to  operator  movement.  More  recently  a  bus  has  been  used 
which  holds  most  of  the  on-shore  optical  and  meteorological 
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equipment  and  employs  the  same  independent  telescope  suspen¬ 
sion  idea.  Ship  motion  originally  made  it  difficult  to  keep 
even  a  diverged  laser  beam  on  the  telescope  long  enough  to 
collect  maningful  MTF  data.  A  gyro-stabilized  tracking  system 
was  designed  to  go  aboard  the  Acania.  At  first  it  was  manned 
by  two  operators  who  sight-aimed  a  gyro-mounted  laser  with 
hand-operated  servo  joy  sticks;  one  controlled  horizontal 
and  the  other  the  vertical  train.  Little  useable  data  was 
obtained  with  this  system  because  of  problems  with  human 
coordination  and  fatigue.  The  solution  was  a  quadrant  detec¬ 
tor  servo  system  which  will  lock  on  a  laser  beam  once  the 
operator  has  aligned  it  with  the  target.  Shipboard  scin¬ 
tillation  detectors  are  also  mounted  on  the  gyro  frame  and 
thereby  maintain  a  constant  relative  position  in  the  incident 
beam  path.  The  gyro  has  been  housed  in  an  aluminum  shed, 
when  aboard  the  Acania,  to  protect  both  equipment  and  operator 
from  the  elements . 

Scintillation  data  was  best  taken  aboard  the  Acania;  a 
detector  with  a  wide  angle  of  acceptance  could  continuously 
sample  a  diverged  beam  from  ashore  more  readily  than  the 
reverse  because  of  the  effect  ship  motion  would  have  on  the 
direction  of  propagation  of  a  beam  originating  on  board. 

The  single-detector  system  originally  employed  required 
excessive  attention  during  data  runs  because  light  from 
environmental  sources  refracts  into  the  beam  path  and  can't 
be  distinguished  from  that  of  the  monitored  laser.  Since  all 
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The  bus  with 
telescope  in 
position.  Transit 
to  the  right  was 
used  in  range-to- 
ship  calculation. 

In  background  stands 
meteorological  tower 
supporting  sensors. 


Naval  Postgraduate 
Research  Vessel  Acania 
with  meteorology  gear 
deployed.  The 
aluminum  shed  which 
shelters  the  gyro 
is  aft  of  the 
pilot  house. 
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Right  side  of 
telescope  showing 
fan^scan  unit, 
laser,  detectors, 
tracking  unit, 
and  associated 
electronic  gear. 


Gyro-stabilized 
tracking  unit.  The 
detectors  are 
mounted  on  top  and 
the  lasers  are 
hidden  from  view 
below. 
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Oscilloscopes , 

Nuclear  Data  Analyzer, 
log-amplifier,  and 
demodulator;  here 
shown  set  up  for 
meteorological  data 
processing  on  the 
bus . 


Dry  lab  on  the 
Acania  with 
scintillation/ 
extinction  data 
processing  gear; 
PIP-400,  PAR 
amplifiers , 
demodulator ,  log- 
converter,  HP9810A 
calculator  and 
plotter.  Oscillo¬ 
scope  &  freq 
analysis  tape 
recorder . 
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photocell  data  passed  through  a  log-amplifier  a  skewed 
distribution  would  result.  In  an  attempt  to  negate  this 
effect  the  beam  had  to  be  interrupted  ten  seconds  out  of 
every  minute  to  zero-out  detector  voltage  generated  by  back¬ 
ground  radiation.  A  dual-detector  system  was  built  with  the 
idea  that  the  output  of  one  pin  diode  detecting  only  background 
could  be  subtracted  continuously  from  that  of  another  xinit 
aimed  at  the  laser  beam. 

Since  the  detectors  (supposedly  identical)  produced 
different  proportional  outputs  extensive  voltage  biasing 
was  necessary.  An  even  greater  problem  was  that  ship  motion 
continuously  changed  the  position  of  the  detector  relative 
to  the  beam  center.  Subsequently  another  single-detector 
system  was  designed,  this  one  employing  a  chopped  laser  source 
and  a  demodulator.  The  chopping  was  accomplished  at  the  source 
by  reflecting  the  laser  light  first  off  of  a  cylindrical 
mirror  which  fanned  the  beam  horizontally  and  then  off  of 
a  flat  mirror  vibrating  at  approximately  two  kilohertz  which 
made  the  beam  scan  vertically.  The  demodulator  unit  took 
the  detector  output  and  electronically  subtracted  the  minima 
(background)  from  the  maxima  (laser  &  background)  and  sent 
the  difference  to  the  log-converter.  Detector  sensitivity 
was  recently  increased  by  replacing  the  pin  diode  with  an 
avalanche  detector.  An  FET  follower  converts  the  high- 
impedance  output  of  the  detector  to  a  low  impedance  which 
can  be  transmitted  a  greater  distance  through  coaxial  cable. 
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This  has  allowed  all  data  processing  gear  on  the  Acania  to 
be  positioned  below-decks  in  the  dry  lab. 

Processing  of  scintillation  data  was  initially  best 
accomplished  through  physical  measurement  of  photographs 
of  resultant  distributions  displayed  on  the  PIP- 400  pulse 
height  analyzer.  Attempts  to  record  and  sxibsequently  analyze 
data  with  integral  methods  on  a  computer  suffered  from  poor 
digitizing.  Digital  output  could  be  obtained  directly  from 
the  PIP-400  by  a  fast  printer  or  from  a  Nuclear  Data  128- 
channel  analyzer  by  teletype.  An  attempt  was  made  to  fit  a 
Gaussian  cxorve  to  the  data  using  least-squares  method.  This 
analysis  proved  highly  successful  when  extraneous  voltages 
did  not  skew  the  experimental  distribution  but  was  painfully 
slow  because  data  had  to  be  manually  input  to  the  computer 
point  by  point.  The  most  recent  improvement  has  been  the 
direct  interfacing  of  the  PIP-400  with  a  compact  computer 
programmed  to  process  scintillation  and  extinction  data  on¬ 
scene  within  six  minutes  of  sampling.  This  was  completed 
by  the  author  and  is  the  svibject  of  this  thesis. 
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III.  THEORY 


A.  SCINTILLATION 

The  theory  of  optical  propagation  in  a  turbulent  atmos¬ 
phere  has  been  extensively  documented  and  those  portions 
employed  by  the  Naval,  Postgraduate  School  optical  propaga¬ 
tion  group  have  not  changed  significantly  since  the  project 
began.  A  thesis  stibmitted  by  B.  C.  Haagensen  [Ref.  3]  in 
June,  1973  contained  an  excellent  outline  of  the  pertinent 
theoretical  derivations  and  an  extensive  bibliography.  The 
author  does  not  intend  to  list  or  verify  familiar  theoretical 
developments  but  rather  to  present  the  significant  end 
products  vital  to  data  processing. 

Theory  states  that  atmospheric  turbulence  is  a  random, 
non-linear,  three-dimensional  phenomenon.  A  turbule  is 
ideally  assumed  to  be  an  isotropic,  homogeneous  pocket  or 
swirl  of  air  with  a  particular  index  of  refraction.  The 
inertial  subrange  delineates  a  region  of  turbule  size  in 
which  refraction  is  theoretically  described;  it  is  bounded 
on  one  side  by  turbules  so  small  (a  few  millimeters)  that 
viscous  dissipation  of  energy  causes  their  immediate  collapse 
and  on  the  other  by  eddies  so  large  (order  of  meters)  that 
they  are  no  longer  homogeneous  and  energy  is  added  to  the 
turbulence  spectrum.  In  the  inertial  siibrange  large  turbules 
form  and  transfer  energy  to  successively  smaller  turbules 
with  minimal  losses  through  viscous  dissipation. 
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Tatar ski  [Ref.  1]  described  the  three-dimensional  power 


spectral  density  of  turbulence,  with  the  equation: 


2  ^-11/3 


(K)  =  0.033  C. 
n 


’N 


where : 


C..  =  refractive  index  structure  constant 


2-ir 

K  =  spatial  wave  number  =  — 
t  =  coherence  length  of  a  turbule 


2 

K  =  tf.  =  lower  limit  of  inertial  subrange 

m  £«  0 


0 


obvioxasly  parameterizes  the  total  amount  of  energy 

in  the  turbulence.  The  random  distribution  of  turbulence 

cannot  be  forecast  and  must  be  described  statistically.  A 

beam  is  modulated  in  a  multiplicative  manner  when  it  passes 

through  successive  turbules  of  differing  index  of  refraction. 

It  follows  that  the  natural  log  of  the  beam  intensity  is 

altered  additively  and  should  follow  a  normal  probability 

distribution  as  a  result  of  the  central  limit  theorem.  C„ 

N 

is  directly  related  to  the  structure  function  of  the  refrac¬ 
tive  index  (the  mean  square  difference  of  the  index  of 
refraction  at  two  points  a  distance  d  apart)  by; 


Dj^(d)  =  <[n(r)  -  n(r+d)]^>  =  d^^^ 
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After  lengthy  derivation  and  a  number  of  assumptions  [Ref. 
4,  pg.  10],  can  be  seen  to  be  directly  dependent  on  the 
variance  of  the  log  intensity  distribution  by; 


=  1.42  z-11/12 


where 


2  2 

=  <[ln(I)  -  ln<I>]  >  =  variance 


=  optical  propagation  path  length 


can  also  be  calculated  by  a  statistical  determination 

2 

of  the  temperature  structure  constant,  ,  since  tempera¬ 
ture  is  the  dominant  factor  in  a  meteorological  formulation 
of  the  index  of  refraction: 


Cjj  =  (79.0  ^  X  10"^)  C, 


where 


P  =  mean  pressure  in  millibars 
T  =  mean  temperature  in  K 
^  2  _  <[T(r)  -  T(r+d)  ]^> 

S  =  ^273 
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A  major  objective  of  the  project  is  to  study  the 
correlation  between  values  obtained  optically  and 
meteorologically . 

B.  GAUSSIAN  MODELLING 

According  to  theory,  periodic  sampling  of  a  laser  beam 
propagating  in  the  inertial  subrange  should  approximate  a 
lognormal  distribution;  if  the  natural  log  of  the  detector 
voltage  (proportional  to  beam  intensity)  defines  abscissa 
values  and  the  ordinate  marks  the  number  of  samples  returning 
a  specific  voltage  the  resultant  plot  should  be  of  the  form: 


2 


where ; 


C  is  a  constant 


X  is  natural  log  of  detector  voltage 


X 


is  the  average  natural  log  of  detector  voltage 


2  2 
a  is  the  variance  (C^j 


Taking  the  natural  log  of  each  side: 


2 


2  ^  V,  ^ 

ax  +  bx  +  c 


In  y  =  - 


X 


2a 


^  ^  t  (InC  - 

a  2a 
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where : 


1 


a 


2 


2a 


b 


X 


2 


a 


c  =  In  C 


2a 


2 


Experimentally  the  natural  log  of  detector  voltage  axis 
is  divided  by  the  PIP-400  into  two  hundred  channels  of  equal 
width  and  the  y  axis  is  the  number  of  counts  (or  times)  the 
equivalent  detector  voltage  (or  beam  intensity)  was  detected 
over  the  sampling  period.  The  number  of  counts  in  each  chan¬ 
nel  after  a  data  run  can  be  assumed  to  be  Poisson-distributed; 
for  one  measurement  the  68%  confidence  interval  is  y  ±  /y. 
Since  the  natural  log  of  y  is  taken  to  obtain  the  quadratic 
form,  the  uncertainty  in  ln(y)  becomes: 


dy  -  -  i 

y  y  ./t 


d[ln(y)  ] 


/y 


So,  for  the  log  distribution,  the  channels  with  the  largest 
numbers  of  counts  have  the  narrowest  confidence  intervals 
(smallest  uncertainty)  and  should  be  weighted  more  heavily. 

The  effect  of  narrowing  confidence  intervals  with  increasing 
counts  is  easily  illustrated  by  plotting  them  on  log  paper. 

In  the  original  program  it  had  been  decided  to  weight  the  data 
inversely  with  the  uncertainty,  i.e.  multiply  by  /y. 
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the:  loendrmrl 
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The  ln(y)  distribution  resulting  experimentally  can  be 
approximated  with  a  polynomial  by  means  of  least-squares 
fitting.  The  method  seeks  to  minimize 
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E  w. [ln(y.)  -  ax.^  -  bx.  -  c]^ 
i=0  ^  ^  ^  ^ 

where: 


=  channel  nvimber  i 
y^  =  counts  in  channel  i 

w^  =  weighting  of  data  in  channel  i,  /y~ 


Setting  the  partial  derivatives  with  respect  to  a,  h,  and 
c  equal  to  zero  and  cancelling  common  terms  gives  three 
equations  with  three  unknowns; 


2  4  3  2 

Ew.x.  Iny.  =  aEw.x.  +  bEw.x.  +  cEw.x. 
11  -^1  11  11  11 


3  2 

Ew.x. Iny.  =  aEw.x.  +  bEw.x.  +  cEw.x. 
ii  i  ii  ii  ii 


Ew.lny.  =  aEw.x.  +  bEw.x.  +  cEw. 

ii  ii  i 


Applying  Cramer's  Rule  [Ref.  6]  to  the  system  determines 
the  coefficients: 
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where 

det (A) 


det  (A^) 


det (A^) 


det  (Aj^) 
^  det  (A) 

det(A2) 
^  det  (A) 

det (A^) 
^  ""  det  (A) 


=  Sw.  X.  ^Sw.x.  ^{Ew.-Zw.x.  } 

3  2  3 

+  Sw.x.  {Zw.x.  Ew.x.-Ew.x.  Zw.} 
11^  11  11111'' 

2  3  2  2 

+  Zw.x.  {Zw.x.  Zw.x.-(Zw.x.  )  } 

11^11  1111  '* 

3 

=  Zw.x.  {  Zw .  X .  Zw .  Iny  .  -Zw .  Zw.x.  Iny  .  } 

11^111^1  1  11  ■^i*' 


2  2  2 
+Zw^x^  {Ew^x^Zw^x^lny^-Ew^x^  Ew^x^  Iny^^} 

2  2 
+  Ew^x^{Ew^x^  Ew^x^lny^-Zw^x^Ew^x^  LhY^} 


=  Ew.x.{Ew.Ew.x. Iny . - Ew . x . Ew . Iny . } 

11  1  11  -^1  111  -'i 


3  2  2 

+  Ew . X .  { Ew . X .  Ew . Iny . - Ew . Ew . x .  Iny . } 


2  2  2 
+  Ew^x^  {Ew^x^Ew^x^  Iny^-Ew^x^  Ew^x^lny^} 
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4  2 

det(A^)  =  Ew.x.  {Lw.x.  Ew . Iny . -Ew. x . Sw. x . 

XX  X  X  xxxx 


3  2  3 

+Ew.x.  {Ew.  X.  Ew.x.  iny..-Zw.x.  Ew.lny.} 

XI  XX  XX  X  11  1  “^1 


+  Ew^x^^  {  Ewj^Xj^^Ew^x^lny^-Ew^Xj^^Ew^x^^lny^ } 


Exactly  the  same  forms  for  the  coefficients  result  from 
applying  Gauss-Jordan  elimination  [Ref.  6].  A  total  of 
eight  distinct  weighted-data  summations  form  the  basic  ele¬ 
ments  which  combine  to  give  the  coefficients.  Coefficient 

"a”  and  a  channel  number-to-log  voltage  conversion  factor 

2 

deliver  the  lognormal  variance  )  which,  with  the  optical 

wavelength  and  range,  determines  the  refractive  index  struc¬ 
ture  constant  •  Coefficients  "a",  "b",  and  the  conversion 

factor  fix  the  average  log  voltage  value  (In 

C.  EXTINCTION 

Extinction  is  the  loss  of  time-averaged  intensity  with 
distance  by  a  beam  propagating  through  the  atmosphere.  No 
extinction  measurements  had  been  taken  by  the  optical  propa¬ 
gation  group  before  September,  1975.  Assuming  an  exponential 
dependence  and  spherical  spreading,  extinction  can  be 
represented  by  an  equation  of  the  form; 

I  =  — 2 —  [-y(R  -  Rq)  ] 


where : 
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I  =  intensity  at  distance  R 


Iq  =  intensity  at  distance  Rq 
Rq  =  some  arbitrary  standard  distance 
R  =  some  arbitrary  distance 
y  =  the  extinction  coefficient 

An  equivalent  form  of  the  equation  is ; 

1  2 

I  =  exp(-yR) 

R 

If  Rq  is  some  small  distance  (on  the  order  of  a  meter) ,  the 
exponential  with  Rq  approaches  one  and  the  above  form  can 
be  simplified  to: 

I  =  -%exp(-yR) 

R 


where 


2 

C  =  IqRq  (a  measure  of  source  intensity) 


Taking  the  natural  log  of  both  sides: 


Inl  =  InC  -  21nR  -  yR 


Solving  for  y : 
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InC  -  21nR  -  Inl 
R 


IR' 


R 


The  constant  C  must  be  determined  at  the  time  of  the  experi¬ 
ment  and  the  time-averaged  intensity  and  range  must  be 
measvired.  Since  detectors  generate  voltages  proportional  to 
the  incident  intensities,  voltage  can  be  substituted  for 
intensity  with  meaningful  results.  Phil  Parrish  investigated 
the  extinction  constant  (C)  and  coefficient  (y)  in  a  con¬ 
current  thesis  [Ref.  10] . 

A  statistical  analysis  was  done  to  determine  the  extinc¬ 
tion  coefficient's  degree  of  sensitivity  to  error  in  C, 

R,  and  I : 


y(C,R,I) 


InC  -  21nR  -  Inl 
R 


iE 

u 


Checking  sensitivity  to  error  with  the  true  value  of  y 

-3 

assumed  to  be  1  x  10  : 
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TRUE  RANGE  form  ^  value 
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Obviously,  errors  in  all  quantities  have  more  effect 
at  shorter  ranges.  Errors  in  the  extinction  coefficient 
appear  to  depend  linearly  on  independent  errors  in  the 
extinction  constant,  time-average  intensity,  and  range  but 
decrease  with  increasing  true  range.  Errors  in  range  have 
a  much  greater  effect  than  equivalent  errors  in  either  of 
the  other  quantities.  The  last  two  calculations  in  the 
previous  table  use  maximxim  error  values  estimated  to  be 
appropriate  with  present  equipment  and  procedure. 
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IV.  INTERFACING 


The  fitting  of  normal  distributions  to  scintillation 
data  through  a  computer  program  using  least-squares  method 
had  proven  highly  successful  but  also  extremely  slow.  The 
PIP-400  pulse  height  analyzer  used  to  accumulate  data  dis¬ 
tributions  had  the  design  capability  to  output  the  information 
through  a  50-pin  plug  in  parallel  BCD  (binary  coded  decimal) 
readout.  The  HP9810A  calculator  in  which  the  processing  took 
place  had  an  optional  parallel  input  interface  card  which 
accepted  BCD  input  through  a  50-pin  plug  at  the  same  voltage 
levels  as  put  out  by  the  PIP-400.  It  was  thought  that  a 
simple  interface  matching  the  proper  lines  would  obviate  the 
pokey  m.anual  input  of  data.  A  number  of  hidden  problems 
were  found  which  made  the  project  much  more  difficult  than 
expected. 

The  first  problem  was  that  both  machines  had  been  designed 
as  masters.  The  analyzer  sent  a  high-to-low  (voltage)  print 
command  and  waited  for  a  low-to-high  flag  from  the  peripheral 
(which  told  it  to  step  to  the  next  channel) .  The  calculator 
sent  a  high-to-low  voltage  command  to  an  external  source 
asking  for  data  and  waited  for  a  low-to-high  flag  which  told 
it  the  information  was  ready.  The  solution  was  to  invert 
both  sets  of  signals  so  that  the  data  request  from  the  calcu¬ 
lator  appeared  as  a  peripheral  flag  to  the  analyzer  and  a 
data  print  command  from  the  PIP-400  appeared  as  an  information 
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ready  flag  to  the  HP9810.  If  the  analyzer  fxanction  knob 
is  set  to  "PRINT"  before  the  calculator  asks  for  data  the 
first  channel  will  be  ignored,  but  since  the  first  (and 
last)  five  channels  of  the  distribution  are  zeroed  by  the 
program  anyway  (because  of  the  high  probability  of  extrane¬ 
ous  coxints  in  these  channels)  the  deletion  has  negligible 
effect.  Since  the  voltage  signal  levels  for  both  machines 
matched,  the  modification  was  completed  with  a  single  hex 
inverter  chip  (Signetics  7404)  using  two  of  the  six  inverters 
available.  A  six  volt  supply  lead  from  the  analyzer  and  a 
100-ohm  resistor  filled  the  inverter's  power  requirement 
(5  ±  .25  volts) . 

Matching  the  data  lines  proved  more  of  a  problem.  Data 
from  the  analyzer  in  parallel  readout  consists  of  six  data 
digits  followed  by  three  address  digits.  Digital  information 
is  relayed  in  BCD  which  requires  four  bitlines  for  each 
decimal  digit.  The  PIP-400  puts  out  a  low  current,  high 
impedance  data  signal  while  the  HP9810  uses  a  pull-up  resis¬ 
tor  to  bring  an  expected  high  current,  low  impedance  signal 
up  to  five  volts  for  a  "1"  bit.  The  calculator  reads  any¬ 
thing  over  2.8  volts  as  a  "1"  bit  and  anything  less  as  a 
"0"  bit.  With  the  command  lines  being  inverted  by  the  chip 
and  the  thirty- four  information  bit  lines  properly  mated 
(inside  a  junction  box  with  two  50-pin  plugs)  a  calculator 
print-out  of  analyzer  data  showed  all  9's;  this  meant  that 
the  pull-up  resistor  in  the  calculator  was  raising  all 
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analyzer  output  signals  to  over  2.8  volts.  This  supposition 
was  substantiated  when  the  correct  data  was  properly  read 
out  of  the  analyzer  by  a  fast  printer  through  the  same  output 
plug. 

The  first  solution  attempted  was  to  install  a  grounded 
quarter-watt  resistor  in  parallel  with  each  signal  line  to 
pull  the  voltage  of  a  "0”  bit  below  2.8  volts.  The  address 
numbers  were  successfully  printed  using  a  5600  ohm  resistor 
on  each  lead  but  the  data  digits  were  decidedly  non-uniform 
in  their  individual  bit  resistor  requirements.  Comparison 
of  calculator  printouts  with  fast  printer  readouts  were  used 
to  determine  the  necessary  resistance;  lead  requirements 
varied  between  4700  and  8200  ohms  for  the  calculator  to  print 
data  digits  accurately.  Variable  (0  -  10,000  ohm)  poten¬ 
tiometers  were  used  to  find  the  range  of  resistance  values 
between  those  which  fixed  a  bit's  output  at  either  "0"  or 
"1".  The  proper  resistance  span  was  smaller,  in  some  cases, 
than  that  separating  successive  standard  resistors  and  was 
found  to  vary  with  time  (possibly  temperature  dependent) . 

Upon  testing  line  voltage  levels  it  was  found  that  while  the 
resistors  brought  a  "0"  bit  from  3.0  to  2.7  volts  it  also 
brought  a  "1"  bit  down  from  5.0  to  2.9  volts  (approximate 
values) .  Since  the  fixed-resistor  interface  was  obviously 
too  finicky  it  was  abandoned. 

The  next  scheme  attempted  employed  operational  amplifiers 
(Signetics  yA747)  constructed  on  silicon  chips  as  unity-gain 
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voltage  followers  and  succeeded  in  generating  the  high 
current,  low  impedance  input  the  calculator  was  designed 
to  use.  Seventeen  op- amps  were  required  to  handle  the  thirty- 
four  bits  of  information  which  made  up  each  data  point.  The 
chips  were  all  connected  in  parallel  with  a  100  mA  dual  power 
supply  (run  on  115  volt  AC)  which  generated  the  ±  15  volts 
direct  current  necessary. 


input  A  . . 


input  B 


■♦1 

14  • 

13* 

•■3 

12  •- 

■♦4 

11  ♦ 

^5 

10 

■♦7 

8 

V  =  -15  Volts 


operational  anplif ier 
socket  wiring  for  tise 
as  a  voltage  follower 


output  A 


(A  &  B  are  BCD  bit  lines) 
output  B 


V  =  +15  Volts 


With  the  op-amps  and  additional  wiring  the  circuitry  was 
too  bulky  to  be  contained  in  the  original  junction  box;  after 
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the  system  had  proven  itself  in  the  field  it  was  redesigned 
and  built  with  a  junction  box  large  enough  to  hold  all  the 
components  internally. 

The  interfacing  of  analyzer  and  calculator  cut  the  time 
to  input  each  data  point  from  seconds  to  milli-seconds; 
with  a  200-channel  analyzer  (PIP-400)  this  meant  a  savings 
of  five  to  ten  minutes  processing  time  on  each  run  and  made 
timely  on-scene  analysis  of  experimental  distributions 
possible . 


36 


The  remnants  of  the  first  interface  using  operational 
amplifiers.  While  successful  it  was  too  vulnerable  in 
an  operational  environment. 


The  interface  in  its  final  form.  To  the  right  are  the  op 
amps  and  inverter  chip,  to  the  left  is  the  power  supply. 
50-pin  plug  on  right  goes  to  calculator,  one  on  left  to 
analyzer;  outlet  plug  is  for  dual  power  supply. 
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V.  COMPUTER  PROGRAM  DEVELOPMENT 


A.  MODIFICATION  OF  THE  ORIGINAL  PROGRAM 

While  mechanically  interfacing  the  pulse  height  analyzer 
(PIP-400)  and  the  computer  (HP9810)  saved  a  great  amount  of 
time  previously  used  to  input  data,  the  computer  program 
which  processed  the  data  was  itself  inefficient.  An  inor¬ 
dinate  amount  of  time  was  spent  labelling,  giving  operator 
instructions,  and  searching  for  subroutines.  The  program 
had  been  written  in  a  straight-line  manner  with  a  minimal 
amount  of  branching  for  recurrent  processes.  The  original 
program  processed  data  for  a  least-squares  fitting  of  a 
Gaussian  and  calculated  the  lognormal  variance  (SIGMA)  and 
the  refractive  index  structiire  constant  (C„)  but  did  not 
compute  the  average  detector  voltage  or  an  extinction  coef¬ 
ficient.  The  program  lacked  the  versatility  to  accept  data 
from  available  analyzers  by  manual  input. 

The  Hewlett-Packard  9810A  calculator  (with  option  003 
installed)  has  a  total  of  2036  locations  in  memory  available 
for  program  steps  of  which  1847  were  used  in  the  original 
program.  Since  the  additional  desired  processing  could  not 
possibly  be  accomplished  in  the  steps  remaining,  a  major 
program  overhaul  was  required.  The  first  measure  taken  was 
to  delete  or  shorten  as  many  labels  and  operator  instructions 
as  possible;  approximately  250  program  steps  were  excised 
in  this  manner  with  no  functional  loss.  More  steps  were 
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saved  by  recording  implemental  constants  on  a  magnetic  card 
and  entering  them  directly  into  the  proper  memory  locations 
immediately  after  the  program  was  loaded. 

The  program  now  had  enough  room  for  additional  processing. 
In  previous  analysis  [Ref.  3,  et.al.)  it  was  assumed  that 
the  natural  log  of  the  average  detector  voltage  (beam  inten¬ 
sity)  was  approximately  equal  to  the  average  natural  log  of 
detector  voltage.  The  validity  of  this  assxamption,  initially 
made  when  computing  the  log- intensity  variance  with  photo¬ 
graphs  ,  is  crucial  in  extinction  measurement  (which  depends 
on  accurate  determination  of  the  average  detector  voltage) . 

The  immediate  concern  was  the  possibility  of  a  functional 
relationship  between  average  detector  voltage  (V)  and  the 
exponential  of  the  average  log  detector  voltage  (exp ( In (V) ) ) 
dependent  on  the  average  value  (ln(V))and  the  standard 
deviation  (a)  of  the  lognormal  distribution. 

A  computer  program  was  developed  which  analyzed  the 
variation  of  the  ratio  of  the  two  quantities  with  increments 
of  either  InCV)  or  0.  The  method  was  to  assume  an  ln(V) 
value  for  a  Gaussian  distribution  and  compute  the  average 
voltage.  The  probability  density  function  of  ln(V)was 
assumed,  therefore,  to  be  of  the  form; 

^  r  (x  -  x) ^ . 
y,n  =  C  exp[ - 5 - ] 

1“  On 


where : 
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X  =  ln(V) 


X  =  In  (V) 


Let  y (V)  be  the  probability  density  function  of  V: 


b  exp(b) 

/  y^  dx  =  /  y(V)  dV 

a  exp  (a) 


dx  =  d{ln(V) }  =  ^ 


^in  =  yin  — 


dV  -^In 


V 


dV 


/  y^Ln^x  =  / 


In(oo)  y 


In 


ln(0) 


V 


dV  =  /  y  (V)  dV 


Since  the  number  of  total  counts  acciamulated  over  a  range 
of  detector  voltage  while  sampling  will  be  identical  whether 
the  display  distribution  is  linear  or  logarithmic: 


y(V)dV  =  Y^^dX  = 


In 


dV 


y(v)  = 


In 

V 


C  „  r  (x-x)  ^ 
=  _  exp  [-  ^  ^ 


2a 


]  = 


C  e^pt.(ln(V)-lniV) 
2a 


V 


Applying  statistics: 

V  =  /vy  (V)  dV 

/y(V)dV 
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For  a  normal  distribution  a  dispersion  index  of  ±3a  gives 
a  confidence  limit  of  99.73%  [Ref.  7,  et.al.)  and  these 
limits  were  imposed  in  a  program  applying  Weddle's  rule 
[Ref.  8]  to  accomplish  the  integration.  The  ratio  of  V 
to  exp(ln(V))  was  found  to  be  independent  of  ln(V)  but  an 
increasing  function  of  a.  With  points  generated  by  the  pro¬ 
gram  (for  a  values  between  0.1  and  1.2  at  increments  of  0.02) 
a  fourth-degree  polynomial  was  found  (using  least-squares 
fitting)  to  generate  average  detector  voltage  values  from 
average  log-intensity.  The  polynomial  coefficients  were 
entered  on  the  magnetic  card  carrying  program  constants 
and  used  to  find  the  average  detector  voltage  after  deter¬ 
mining  InCV)  (from  the  coefficients  of  the  least-squares 
polynomial  fitting  data  to  a  lognormal  distribution) .  An 
easier  method  was  subsequently  employed  which  can  be 
justified  by  proceeding  further  with  the  preceding  derivation. 

^  ^  /Vy(V)dV 
Jy  (V)  dV 


but 


V  =  exp(ln(V))  =  exp(x) 


and 


y(V)dV  =  y^^  dx 
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fC  exp(x)  exp [ 


V  = 


/exp(x) 

/yina>= 


(x-it)  ^ 

2a^ 


]  dx 


/C  exp  — ]  dx 


2a 


After  combining  exponentials  and  completing  the  square  in 

_  2 

the  nximberator,  and  defining  x'  =  x  -  a  : 


C  exp{x+;^)  /exp  ^ — ]  dx 


(x-x') 


V  = 


2a 


( 

C  /exp  A- — ]  dx 

2a"^ 


The  integrands  in  numerator  and  denominator  are  Gaussian 
distributions  with  the  same  standard  deviations  and, 
therefore,  the  same  areas. 


V  =  exp{x  +  =  exp  ( In  (V))  exp(^) 


This  derived  result  can  be  found  in  meticulous  statistical 
literature  [Ref.  9].  Since  ln(V)  and  a  were  already  avail¬ 
able  in  the  program  the  necessary  calculation  was  easily 
included  and  required  significantly  less  memory  than  the 
fourth-degree  polynomial  conversion.  Incidentally,  it  was 
found  that  the  polynomial  method  was  accurate  within  2% 
for  standard  deviation  values  between  0 . 1  and  0.9! 

With  an  assxamed  extinction  constant  the  calculation  of 
an  extinction  coefficient  is  easily  made,  i.e.; 

In(^) 
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where : 


U  =  extinction  coefficient 

2 

C  =  extinction  constant  (V^R  ) 

o  o 

V  =  average  detector  voltage 

R  =  range 

This  computation  was  added  to  the  program  with  the  extinction 
constant  becoming  one  of  the  operator  input  parameters. 

The  next  program  modification  was  the  addition  of  a 
section  which  could  process  linear  data  to  find  the  average 
detector  voltage  and  extinction  coefficient.  The  main 
purpose  was  to  correlate  results  with  those  from  the  log¬ 
arithmic  section  but  this  was  not  done  to  any  great  extent. 
The  average  channel  number  is  determined  by  dividing  the 
summation  of  channel  number  times  channel  counts  by  the 
summation  of  counts.  The  average  channel  number  is  then 
modified  in  a  linear  scaling  process  to  give  the  average 
voltage.  The  extinction  coefficient  was  computed  in  exactly 
the  same  way  as  in  the  logarithmic  section. 

After  the  linear  portion  was  added  the  program  almost 
filled  available  memory  and  was  found  to  execute  more  slowly 
than  before  —  especially  plotting  data.  The  program  showed 
minimal  utilization  of  subroutines  and,  consequently,  the 
program  wasted  memory.  The  subroutines  used  were  scattered 
throughout  the  body  of  the  program  which  meant  wasted  time 
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spent  searching  after  branching  statements.  The  computer 
still  took  data  input  only  from  the  PIP-400  through  the 
interface.  A  Nuclear  Data  128  channel  analyzer  which  inter¬ 
faced  with  a  teletype  was  available  as  a  back-up  but  the 
program  had  to  be  modified  to  accept  manual  input  and  a 
different  channel  capacity.  With  all  the  discrepancies 
a  major  reorganization  was  needed. 

The  HP9810A  calculator  branches  to  subroutines  by  going 
to  the  program  beginning  and  searching  through  until  finding 
the  called  label.  This  requires  an  inordinate  amount  of 
time  to  complete  ”do-loops"  and  subroutines.  This  problem 
was  solved  by  making  processes  which  were  reiterative  or 
used  by  more  than  one  section  of  the  program  into  subroutines 
and  placing  them  in  front  of  the  main  body.  "Flag”  state¬ 
ments  were  used  to  allow  automatic  or  manual  data  input  from 
the  PIP-400  or  Nuclear  Data  analyzer  respectively. 

The  original  program  used  a  large  amount  of  memory  in  a 
calculation  to  filter  out  bias  introduced  in  the  scintillo¬ 
meter.  When  the  single-detector  system  for  a  chopped  laser 
source  was  first  put  together  the  demodulator  would  output 
a  few  tenths  of  a  volt  with  no  input  present  thus  biasing 
the  sample  distribution.  This  offset  had  to  be  taken  into 
account  when  processing  data  because  of  its  effect  when 
using  the  logarithmic  scale.  In  subsequent  demodulators 
the  offset  problem  was  solved  and  the  calculation  became 
extraneous.  When  the  process  was  removed  from  the  program 
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the  only  discernable  difference  was  an  increased  speed  in 
plotting  data. 

The  net  result  of  all  the  program  modifications  pre¬ 
viously  discussed  was  a  drop  in  the  average  processing 
time  of  an  experimental  distribution  from  over  ten  minutes 
to  less  than  six  minutes. 

B.  METEOROLOGY  PROGRAMS 

Just  prior  to  the  most  recent  R/V  Acania  experiments 
(May,  1976)  Professor  Schacher  (associated  with  the  atmos¬ 
pheric  turbulence  group)  requested  a  computer  program  to  fit 
a  Gaussian  curve  to  data  manually  input  from  the  Nuclear 
Data  pulse  height  analyzer.  No  explanation  of  the  relevant 
theory  was  given  or  required  to  distill  the  necessary 
processing  from  the  existing  scintillation/extinction  program. 

Upon  analyzing  the  processed  data  it  was  found  that  the 
Gaussian  fit  was  a  poor  one.  Prof.  Schacher  suggested  the 
resultant  data  distributions  appeared  more  Lorentzian  and 
requested  another  program  fitting  data  to  that  distribution. 

The  Lorentzian  distribution  is  of  the  form; 


A 


2  —  2 
r  +  (x-x) 


2  —  —2  2 
X  -  2xx  +  X  +  r 


1 

y 


This  can  again  be  fitted  by  least-squares  method  to  a 
quadratic  and  the  coefficients  evaluated  to  get  the  average 
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voltage  (after  scaling  x,  the  average  channel  number)  and  r 
(after  scaling) .  The  weighting  scheme  differed  from  that 
used  for  the  Gaussian.  The  uncertainty  is  (again  assuming 
the  nximber  of  counts  in  each  channel  is  Poisson-distributed) 

,,lv  -2,  -dy  -t/y  -3/2 

d(-)  =  -  y  dy  =  — f  =  -y 

^  y  y 

Again  the  data  was  weighted  inversely  with  its  uncertainty 
(multiplied  by  the  number  of  counts  in  the  channel  to  the 
3/2  power) .  The  least-squares  fit  attempts  to  minimize: 

Zw.(-^  -  ax.^  -  bx.  -  c) 

1  y.  1  1 


where 


bJ 


i 


3/2 


Y .  =  number  of  counts  in  channel  i 
1 


X.  =  channel  niomber  i 
1 


a 


1 

A 


b 


-2x 

A 


C 


+ 

A 


r 


2 


Since  all  the  machinery  of  the  least-squares  fitting  of 
a  quadratic  applying  Cramer's  rule  already  existed,  it  was 
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a  relatively  simple  matter  to  modify  the  existing  Gaussian 
program.  The  comparative  results  definitely  point  to  a 
Lorentzian  distribution,  which  was  not  predicted  by  theory. 
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ENTER  CALI  BRA' 
CALIBRATION  V( 

TION  CHANNELS, 
OLTAGES 

CONSTRUCT  SCALE  FACTOR 

ENTER  EXTINCTION  CONSTANT^ 
WAVELENGTH,  RANGE,  PAR  GAIN, 
MAXIMUM  COUNT,  TIME 

ENTER  DATA  AND  PLOT. 

CON- 

STRUCT  SUMMATIONS  FOR 

LEAST- 

SQUARES  FIT. 

CONSTRUCT  COEFFICIENTS  OF 
LEAST-SQUARES  QUADRATIC 


CONSTRUCT  LOGARITHMIC 
VOLTAGE  variance 


PLOTTER  PRINT;  DATE,  TIME, 
I/1BD  (wavelength),  RNGE, 
(range),  $  SQ  (variance), 

$  (standard  deviation) 


construct  and  plotter  print 
CNSQ  (c^  squared),  CN  (C|^) 


construct  and  plotter  print 
AVE  V  (average  voltage) 


A  SIMPLIFIED  FLOW  CHART 
OF  THE  NEW  PROGRAM 


YES 


CALIBRATION  POINTS 
PRINTED  OUT 


LINEAR 


ENTER  DATA,  CONSTRUCT  SUMMATIONS 
FOR  AVERAGE  VOLTAGE 
-  -'1  '  ' 


CONSTRUCT  AND  PRINT  AVE  V 

(average  voltage) 


CONSTRUCT  AND  PRINT  EXT  MU 

(extinction  coefficient) 


I  END^ 


CONSTRUCT  AND  PLOTTER  PRINT 
EXT  MU  (extinction  COEFFICIENT) 
PLOTTER  PRINT  GAIN  (PAR  GAIN) 


PLOTTER  PLOT  COMPUTED  GAUSSIAN 
CURVE  AND  CALIBRATION  MARKERS 


I  END^ 
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001034  *000000  *000001  *000000 
000000  *000000  *000000  *000000 
000000  *000000  *000000  *000000 
000000  *000000  *000000  *000000 
006291  *008860  *005971  *000000 
000000  *000000  *000000  *000000 
000000  *000000  *000000  *000001 
000000  *000000  *000000  *000000 
000000  *000000  *000000  *000000 
006638  *003506  *000000  *000000 
000000  *000000  *000000  *000000 
000000  *000000  *000000  *006799 
000000  *000000  *000000  *000000 
000000  *000000  *000000.  *000000 
000000  *000000  *000000  *000000 
000000  *000000  *000000  *000000 


*000000  *000001  *000000  *000000 
*000000  *000000  *000000  *000000 
*000000  *000000  *000000  *000000 
*000000  *000000  *000000  *000000 
*000000  *000000  *000000  *000000 
*000000  *000000  *000000  *000001 
*009367,  *016270  *000002  *000000 
*000000  *000000  *000000  *000000 
*000000  *000000  *000000  *000007 
*000000  *000000  *000000  *000000 
*000000  *000000  *000000  *000000 
*002922  *000000  *000000  *000000 
*000000  *000000  *000000  *000000 
*000000  *000121  *007691  *003369 
*000000  *000000  *000000  *000000 
*000000  *000000  *000000  *000000 


An  example  of  the  teletype  output  of  the  voltage  calibration 
channel  numbers  used  in  scaling  (reproeduced  from  1  March  1976) . 
Notice  the  extraneous  counts  collected  in  channel  1.  The 
points  would  be  entered  as  follows; 

Voltage  Channel  #  Change  in  # 


0.25 

32.97 

0.50 

52.63 

19.60 

1.00 

72.35 

19.72 

2.00 

91.30 

18.95 

4.00 

110.30 

19.00 

PICTURE  OF  EXPERIMENTAL  DISTRIBUTION  AS  DISPLAYED  ON  THE 
PIP-400  PULSE  HEIGHT  ANALYZER 
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B 

B 


oinihijiwjnw-H 
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EXAMPLE  OF  PLOTTER  OUTPUT  FROM  THE  OLD  PROGRAM 


i 
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EXAMPLE  OF  PLOTTER  OUTPUT  FROM  THE  NEW  PROGRAM 


I 


SET  FLAG  FOR 

MONTH? 

CAL  RHIWrLUT 

4. 

11701.073 

DAY? 

3117.074 

30. 

7040.103 

YEAR? 

10768.104 

76. 

IO883.134 

FLAG  -  CAL  PTS 

10939.164 

1ST  CAL? 

8756.194 

61.0 

2194.195 

2ND  GAL? 

EWTSR  IST  CAL 

117.0 

EWTER.  2Na  CAL 

3RO  CAL? 

EliTER  3RD  CAL 

173.0 

ENTER  VI  IN  Y 

VI -Y  AND  V3-X 

AND  V3  IN  X. 

0.50 

ENTER.  VOLTAGE 

8.00 

OFFSET- 

EXTINCT  CONST? 

ENTER.  WAVELENGTH 

15300.00 

IN  MICRONS 

WAVELENGTH 

ENTER-  RANGE 

(MICRONS)? 

IN  METERH 

0.6328 

ENTER-  DATE 

RANGE  (METERS)? 

ENTER  TIME 

1625. 

SET.  ANALYZER-  TO 

GAIN? 

PRINT.  OUT.  AND- 

200. 

SE'IL  FLAG  FOR 

MAX  COUNT? 

DATA  PRINT  OUT 

8000. 

ENTER  MAX  COUNTL 

TIME? 

1311. 

FLAG  -  LINEAR 
FLAG  -  LOAD  DATA 

DATE 

PIP  ON  PRINT? 

20776. 

TIME 

1245. 

END. 

kJAVELENGTH 

MONTH? 

0.6328 

4. 

RANGE 

DAY? 

900. 

28. 

YEAR? 

183.092 

76. 

291.093 

1ST  CAL? 

503.094 

2NDl  gal? 

569.095 

3RD  GAL? 

645.096 

VI -Y  AND  V3-X 

640.097 

MAX  COUNT? 

579.098 

80000.00 

758.099 

TIME? 

714.100 

1555.00 

767.101 

LOAD  DATA 

849.102 

(etc) 

END 

Reproductions  of  computer  printer  tapes  made  during  data 
processing:  original  program  with  cal  and  data  printout 
(left),  present  program  (upper  rt.),  meteorology  program 
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s 

El 


El 

m 

El 

ID 


El 

J 


-  oi  m  K  ID  M  J  ri  w  Ei 
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EXAMPLE  OF  PLOTTER  OUTPUT  FROM  THE  METEOROLOGY  PROGRAM 


s 

s 
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EXAMPLE  OF  PLOTTER  OUTPUT  FROM  THE  METEOROLOGY  PROGRAM 


APPENDIX  A 


The  following  lists  the  connections  made  in  mating  the 
output  plug  of  the  PIP-400  to  the  input  plug  of  the  interface 
for  the  HP9810A; 

CALCULATOR 

PIP  LEAD  #  LEAD  #  COLOR 


LEAD  FUNCTION 

DATA  c 

10^  -  bit  1 
"  2 
"  4 

4  " 

10^  -  1  bit 

2  " 

4  " 

8  " 

ADDRESS  - 

10  -  1  bit 

2  " 

10^  -  1  bit 
2  " 

4  '• 

10^  -  1  bit 
2  " 

4  " 

8  " 

DATA  ^ 

10^  -  1  bit 
2  " 

4  " 

2  8  " 
10  -  1  bit 

2  "  ■ 
4  " 

1  3  " 

10-^  -  1  bit 
2  " 

4  " 

8  ** 

10°  -  1  bit 
2  " 

4  " 

8  " 

Print  Command 

Hold  Command 

Ground 

+  6  volts 


1 

6 

2 

7 

3 

8 

4 

9 

5 

10 

6 

11 

7 

12 

8 

13 

9-12 

N/C 

13 

30 

14 

31 

15,16 

N/C 

17 

34 

18 

35 

19 

36 

20 

37 

21 

38 

22 

39 

23 

40 

24 

41 

25 

14 

26 

15 

27 

16 

28 

17 

29 

18 

30 

19 

31 

20 

32 

21 

33 

22 

34 

23 

35 

24 

36 

25 

37 

26 

38 

27 

39 

28 

40 

29 

red 

white/black 

white/green 

blue 

white/blue 

green 

white/gray 

yellow 


white 

white 

white 

II 

II 

It 

white 

II 

II 

II 

white/ violet 
white/black ( 2 ) 
white/yellow/black 
white/blue/black 
white/red/black 
white/orange/black 
black 

white/gray /black 

white/yellow 
white/brown/black 
white/orange 
white/green/black 
white/ violet/black 
white/red 
white/brown 
white/black 
white 


41  inverted  49 

42  "  47 

43  chassis  and  power  supply  white 

44-48  N/C 

49  hex  inverter  white 

50  .  N/C 


N/C  designates  leads  which  are  not  connected 
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APPENDIX  B 


The  following  is  a  listing  of  the  data  storage  registers 
utilized  by  the  program: 


REGISTER 


USE 


000 

001 

002 

003 

004 

005 

006 

007 

008 

009** 

010 

Oil 

012 

013 

014 

021* 

022* 

023* 

025 

028* 

030 

031* 

041* 

050 

051* 

052* 

053* 

054 

055* 

056* 

057 

058 

060 

061 

065 

068 

069 

078* 

079* 

095 

097 

098 

099** 

100** 

102** 

103 

104 


computed  average  logarithm  of  detector  voltage 


computed  summations  for  least-squares  fit;  001, 
002  and  003  are  used  subsequently  to  store  the 
coefficients  of  the  least-squares  polynomial 


9999 

determinants  (formed  from  s\immations)  which  give 
coefficients  of  the  least-squares  polynomial 

computed  standard  deviation  of  least-squares 

Gaussian  fitted  to  data 

first  calibration  channel  number 

second 

third 

computed  variance  of  least-squares  Gaussian, 

later  used  again  in  Cj^  calculation 

time 

wavelength  in  meters 

range  in  meters 

wavelength  in  microns 

used  in  calculation 

month 

day 

year 

5/max  count* 
gain 

extinction  constant 
computed  average  detector  voltage 
computed  extinction  constant 
data  point 

data  point  weighting 

channel  number-to-voltage  scaling  constant 
do-loop  accumulators  for  the  linear  calculation 
of  average  detector  voltage 
first  calibration  voltage 
third 

varying  pen  coordinate  value  for  axis  labelling 
incremental  pen  coordinace  values  for  axis  labelling 
label  values  in  axis  labelling 
236  (range  of  x-axis) 

11.8  (range  of  y-axis) 
alpha  scale  factor  (0.1,  FMT  11) 
computer-scaled  x-coordinate 
"  y-coordinate 
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REGISTER 


USE 


105** 

106** 

107** 

108** 

a  &  b 


10*9  (followed  by  FMT  13) 

—  1  fl 

218  (followed  by  FMT  12) 

'  max^ 

used  many  times  through  out  the  program 


*  marks  those  registers  whose  values  are  input  by  the  operator 
**  marks  those  registers  with  values  stored  on  a  magnetic  card 


The  following  is  a  list  of  label  titles  in  the  order 
they  appear  in  the  program: 


LABEL 

A* (UP,DN) 

B*(l,2) 

C*(3,4) 

D*(A,B,C,5) 

E*  (F) 

F* 

G* 

H*(6,7) 

I* 

J* 

K* 

L(K) 

M(N,H, J,D, 
G,I,E) 


N(J,A,B,I, 

M,8,9) 


FUNCTION 

takes  manual  or  interface  input;  puts  data  in 
060,  channel  in  a,  5/max  count  in  b 
puts  channel  niomber  in  y  and  upper  data  channel 
limit  in  x 

puts  channel  number  in  y  and  upper  channel 
limit  in  x 

plots  data  points  and  constructs  summations 

necessary  for  least-squares  fitting 

plots  the  computed  Gaussian  and  calibration  points 

plots  marker  at  calibration  point 

constructs  a  special  floating  point  form  of 

a  nxamber  to  be  plotted 

labels  plotted  x  &  y  axes 

computes  the  extinction  coefficient 

prints  operator  instructions : 

FLAG  -  LOAD  DATA;  PIP  ON  PRINT? 
lists  calibration  points 

preliminary  program  set-up  (month,  day,  year, 
calibration  points,  extinction  constant) 
completes  program  set-up  (wavelength,  range, 
gain,  max  count,  time);  branches  to  label  N 
for  linear;  plots  and  labels  x  &  y  axes;  plots 
data;  calculates  desired  constants;  plots  con¬ 
stants,  computed  Gaussian,  and  calibration 
points;  returns  to  label  M 

linear  portion  of  program;  computes  average 
detector  voltage  and  extinction  constant, 
and  prints  them,  returns  to  label  M 


*  marks  labels  which  are  subroutines 
0  hold  subroutines  and  labels  branched  to 
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APPENDIX  C 


The  following  is  a  complete  listing  of  the  program  as 
printed  out  by  the  calculator; 


0000 

GTO 

44 

0050 

LBL 

51 

0100 

+ 

33 

0001 

LBL 

51 

0051 

DN 

25 

0101 

DN 

25 

0002 

L 

72 

0052 

LBL 

51 

0102 

CNT 

47 

0003 

LBL 

51 

0053 

UP 

27 

0103 

S/R 

77 

004 

A 

62 

0054 

FMT 

42 

0104 

LBL 

51 

0005 

IFG 

43 

0055 

1 

01 

0105 

B 

66 

0006 

SFL 

54 

0056 

8 

10 

0106 

a 

13 

0007 

GTO 

44 

0057 

0 

00 

0107 

UP 

27 

0008 

LBL  . 

51 

0058 

1 

01 

0108 

IFG 

43 

0009 

UP 

27 

0059 

XTO 

23 

0109 

SFL 

54 

0010 

FMT 

42 

0060 

+ 

33 

0110 

GTO 

44 

0011 

1 

01 

0061 

0 

13 

0111 

LBL 

51 

0012 

8 

10 

0062 

0 

13 

0112 

1 

01 

0013 

3 

03 

0063 

UP 

27 

0113 

CNT 

47 

0014 

FMT 

42 

0064 

CLX 

37 

0114 

1 

01 

0015 

3 

03 

0065 

STP 

41 

0115 

9 

11 

0016 

3 

03 

0066 

XTO 

23 

0116 

5 

05 

0017 

21 

0067 

6 

06 

0117 

GTO 

44 

0018 

UP 

27 

0068 

0 

00 

0118 

LBL 

51 

0019 

CLX 

37 

0069 

FMT 

42 

0119 

2 

02 

0020 

X=Y 

50 

0070 

1 

01 

0120 

LBL 

51 

0021 

GTO 

44 

0071 

8 

10 

0121 

1 

01 

0022 

LBL 

51 

0072 

3 

03 

0122 

CNT 

47 

0023 

A 

62 

0073 

LBL 

51 

0123 

1 

01 

0024 

CNT 

47 

0074 

DN 

25 

0124 

2 

02 

0025 

EEX 

26 

0075 

XFR 

67 

0125 

3 

03 

0026 

CHS 

32 

0076 

5 

05 

0126 

LBL 

51 

0027 

3 

03 

0077 

4 

04 

0127 

2 

02 

0028 

X 

36 

0078 

XTO 

23 

0128 

S/R 

77 

0029 

DN 

25 

0079 

b 

14 

0129 

LBL 

51 

0030 

UP 

27 

0080 

1/X 

17 

0130 

C 

61 

0031 

INT 

64 

0081 

UP 

27 

0131 
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